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Abstract 

We investigate the role of unitarization effects in virtual photon-proton 
(7*p) interactions at small x. The (/(/-fluctuation of the initial photon is sepa- 
rated into a small distance and a large distance component and a model for 
the unitarization of each component is proposed. The Born approximation 
for the small size component is calculated using QCD perturbation theory. 
Reggeon diagram technique is used in order to obtain a self-consistent scheme 
for both total "f*p cross section and diffractive production. The model gives a 
good description of HERA data in the small-x region, with a single Pomeron 
of intercept 1.2. 
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1 Introduction 



The present work is an extension of our previous one [|TJ on the investigation of 
unitarity effects in small- x processes. It was found experimentally at HERA that 
both the total cross section of a highly virtual photon, a^*p , and the cross section 
for its diffractive dissociation have a fast increase with energy. This is related to a 
fast increase of densities of quarks and gluons as the Bjorken variable x decreases. 
The dynamics of such very dense partonic systems is very interesting and has been 
studied by many authors both in deep inelastic scattering (see ref. for reviews 
and ref. |3) for some recent papers) and in high energy nuclear interactions 
Unitarity effects should stop the increase of densities at extremely small x and lead 
to a "saturation" of parton densities. It is important to determine the region of x 
and Q 2 where the effects of saturation become important. 

We study this problem using reggeon calculus M with a supercritical Pomeron 
(ap(0) — 1 = A > 0) and the partonic picture of 7*p-interaction in QCD. In our 
previous paper |] we used this approach for the description of HERA data in the 
region < Q 2 < Ql (Q 2 , ~ 10 GeV 2 ), where the effects of unitarity are most impor- 
tant. It was shown that, with a single Pomeron of intercept 1.2 and multipomeron 
exchanges (unitarity effects), it is possible to obtain a self-consistent, simultaneous, 
description of both the total 7*p-cross section and diffractive production in high- 
energy 7*p-interactions. In such approach, it is convenient to consider the process of 
7*p-interaction in the laboratory frame as an interaction of the gg-pair, produced by 
the photon, with the proton. We separated the gg-pair fluctuation into two compo- 
nents, "aligned" component with a strongly asymmetric sharing of the momentum 
fraction z between q and q, and the rest ("symmetric") component. Such a sepa- 
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ration is important at large Q 2 , where the first component has a large transverse 
size, while the "symmetric" component has a size r ~ 4 and thus has a small cross 
section ~ 1/Q 2 of interaction with the target. Both components give a contribution 
to the a!y?p which behaves as 1/Q 2 at large Q 2 , but the "aligned" component gives 
the main contribution to the diffraction production cross section. Triple Pomeron 
diagrams were also included in our model. 

In this paper we propose a more direct separation of the two components of the 
gg-pair, which is valid also for small Q 2 . The separation into a small size (S) and 
a large size (L) components of the qq pair is now made in terms of the transverse 
distance r between q and q. The border value, tq, is treated as a free parameter - 
which turn out to be r ~ 0.2 fmQ. 

For the S"-component, with r < r , we use the expression for the 7*p total cross- 
section obtained in perturbative QCD MM. 



a { X )T(L \s,Q 2 ) = j Q d 2 r J o dz ^ (L) (r,z,Q) " a 5 (r,s,Q 2 ) , (1) 

where T and L correspond to transverse and longitudinal polarizations of a virtual 
photon, ip T ^ (r, z) are the corresponding wave functions of the gg-pair: 



ro „ r 1 
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V(r,z,Q) 2 = ^?£^ 2 + (i _ z 2 )]e 2 Kl{er) +m 2 q K 2 (er)}, (2) 



An 2 „ 



and 



TP L (r,z,Q) 2 = Q -^T / e 2 q {AQ 2 z 2 (l-z) 2 K 2 (er)}, (3) 

where e 2 = z(l — z)Q 2 + m 2 . i^" an d ^i are McDonald functions. The sums are 

over quark flavors and we have taken m u = rrid = m s = ms- 

^This value agrees with the correlation length of nonperturbative interactions observed in lattice 
calculations jql. 
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<?s(r, s, Q 2 ) is the total cross section for the interaction of the qq-p&ir with the proton. 
For the interaction of a small size dipole 

a s (r,s,Q 2 )=r 2 f(s,Q 2 ) . (4) 

As for the L component, we use the same parametrizations introduced in ref. jl| for 
the aligned component. 

2 The model 



We write the 7*p total cross section 



4f^Q 2 ) = ^^ F ^Q 2 ) , (5) 

in the following form, using the impact parameter (b) representation 



<r%\s, Q 2 ) = 4 j d 2 b a { X\b, s, Q 2 ) (6) 



4f(b, s, Q 2 ) = g 2 L (Q 2 ) atHb, s, Q 2 ) + 4 tot) (6, a, Q 2 ) . (7) 

The function g\{Q 2 ) determines the coupling of the photon to the large size qq 
pair and is chosen in the form ffl 



9l(Q 2 ) = ^% (8) 

where g1(0) and m 2 L are phenomenological parameters. 

The cross section for the L-component, cr^ ot \ in the impact parameter space, is 
chosen in the quasi-eikonal form || 



ar\b,s,Q^= 1 ' eM ' C ^ S ^ Q2)) , (9) 



x ^ Q ' ) = TTfWw +xUb ^ ■ (10) 

The eikonal functions Xlo = ^ f) are written in a standard Regge form 



where 
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A fc = a*(0)-1 , t = £n^-Q- , = + aU . (12) 

so + v 

Here «fc(0) is the intercept of trajectory k and a' k its slope. The values of the 
radii Rl kL , based on ref. |TIJ, are given in Table 1. The quantity £ is chosen in such 



a way as to behave as £n^ for large Q 2 and as £n-^ for Q 2 = 0. 

The coefficients C[ and C*/ determine respectively the residues of the Pomeron 
and /-reggeon exchanges in the gg-proton interaction. The coefficient C — 1.5 takes 
into account the dissociation of a proton ||. 

We turn next to the denominator of eq. fllOl). The constant a si given by 
( ._ , where 5^(0) is the proton-Pomeron coupling and rppp(0) is the 



_ _ 3p P (°) rppp(0) 



triple Pomeron coupling, both at t — 0. The function X3{b, s,Q 2 ) is given by eq. 
(|3lD of section 3. 

With a = 0, the model described above is a standard quasi-eikonal model with 
Born terms given by Pomeron plus / exchanges. The denominator in eq. (^) corres- 
ponds to a resummation of triple Pomeron branchings (the so-called fan diagrams). 
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(For a full discussion on the interpretation of this denominator see ref. [|TJ). Thus, 
expressions (|9]) and (|I0) correspond to a sum of diagrams of the type shown in Fig. 1. 

We turn next to the S component. In this case we put, in complete analogy with 
eqs. (D-flllD 



ar\rAs^)= l - eW{ ~ C X " {rASm) , (13) 



w(rA.,<n = T ^*£> , (14) 
1 + a X3{b,s } Q 2 ) 



te » (r ' 6 '«=S5 exp ( Apf -«£«) • (15) 

with A^ P = R 2 PS + a' P £. 

Note that the contribution of the /-exchange to the S component is very small 
and has been neglected 0]. The condition (f|), valid for fixed s and Q 2 as r — > 0, is 
a property of the single Pomeron exchange. Thus a factor r 2 has been introduced 
in eq. (|T5|). 

Finally os(r, s, Q 2 ) in eq. (P is obtained from o"5-(r, b, s, Q 2 ), defined by eqs. (|T3|) 
to ([T5|) , as (see eq. (0)) 

a 5 (r, S ,Q 2 )=4|rf 2 6a 5 (r,6, S ,Q 2 ) . (16) 

Inserting this expression in eq. (|]) we obtain the transverse and longitudinal con- 
tributions of the ^-component to the total j*p cross-section. 
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3 Diffractive production 

Following ref. [JXJ] we express the total diffractive dissociation cross-section of a 
virtual photon as a sum of three terms 



a 



(diff) 
7*P 



i=L,S 



where 



(17) 



aP=4gl(Q 2 ) (a^\b,s,Q 2 )Yd 2 b 



J Jo Jo 



^ L (z,r)\ 2 (af(r,b,s,Q 2 )) 2 (19) 



a PPP = 2g 2 L (Q 2 ) J x L pp P (b,s,Q*) C -«^(M,0») d % 

+2 / d 2 b f l d 2 r f 1 rfzy|f L (v) 2 Xpp P (b,s,Q 2 ) e - 2 ^(^,Q 2 ). (20) 

J Jo Jo 7Z~t 1 



Here 



X L pp P (b, s, Q 2 ) = a X P L (b, s, Q 2 ) X s(b, s, Q 2 ) (21) 

and 

Xp PP (r, b, s, Q 2 ) = a Xs (r, b, s, Q 2 ) X3 (b, s, Q 2 ) (22) 

where XL(b, s,Q 2 ) is given by the first term of eq. (|10D and X3(b,s,Q 2 ) is defined 
by eq. (|31|). Using this expression, we see that, to first order in a, <Jppp consists of 
the sum of a triple Pomeron (PPP) term plus a PfP one. We call this sum triple 
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Pomeron, although the second one is an interference term. For the total diffractive 
production cross-section, that includes the diffraction dissociation of a proton, eqs. 
(|l8|)-(p0|) must be multiplied by the same factor C—1.5 of the total j*p cross-section. 

At HERA, differential diffractive cross sections are given as a function of /3 — 
m ? + q2 , where M is the mass of the diffractively produced system, or of xp = x//3. 

(3) 

They are usually integrated over t, and the function F^jj is introduced 



P 2D An 2 a e , m , I P dxpdt ^ ^ 

In our model, this function can be written as a sum of three terms 



Here 



F?J = ( £ mi*, Q 2 ' P) + F^pppix, Q\ (3) ] . (24) 

i=L,S 



,(3) Q 2 9l(Q 2 ) V- f,2,..i.. k ^ i+ A fe -A /(l _^n H Q 2 ) 



xpF 2DL - L J J2 J ^Xlxl f3 max dpn Ai+Ak _ A ( , _ B)npm ( 25 ) 



and 



P 2DS 47ra e . m .\ s f CZ - 2/3)2 s fffa f ^(l - f3) ) 

where (3 = ffijfr, (3 m in = ^Mss = 10x and f3 max = M ^ +Q 2 with M 2 min = 4m£. In 
eq. fl2"5[) <jf^ B corresponds to eq. (|18|) keeping only the linear term in af ot ^ and 
Xl 1S the contribution of the P(f) in eq. (|10|). 

The /3-dependence of the S'-component has been taken from the QCD results of 
ref. [p| . The /3-dependence of the L-component was chosen according to ref. [ 12 
and 
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B 'C« = -5 + i(^] • (27) 



with c = 3.5 GeV 2 . 

The triple-Pomeron (i.e. PPP plus PfP) contribution, F2p,ppp(x,Q 2 , P), is 
given by 

xp F 2 % PP (x, Q 2 , (5) = x P F^ PP (x, Q 2 , (3) , (28) 

°ppp 

where a PPP is given by eq. (|20"D, its Born term, (J PPP , by the same equation with 
C — 0, and 



Here 



x p F 2Dppp(x,Q 2 ,(3) 



Q' 2 



An 2 a P 



2a d 2 b X3 (b,s,Q 2 ,P)x 



rro rl o 

^(Q 2 )xf(M,Q 2 ) + E/ *r dz^ L (r,z) X s(r,b,s,Q 

JO JO 



(29) 



X 3 (s,b J Q 2 J (3) = 7fcexp 
fe=p,/ 



6 2 I fp\ Ak (1 - /5fHQ 2 )+4 



V 



4A fc (| 



(30) 



where 7p = 1, 7/ determines the strength of the P/P-contribution relative to the 
PPP one, and = R\ k + a' k £n Ogj . The function X3( s > b, Q 2 ), which enters in eqs. 
(ID, O, (P) and (H) is given by 



x 3 M,Q 2 ) = / -^x 3 M,Q 2 ,/?) • 



(31) 



Since the triple Pomeron formula is not valid for low masses, we use here M„ 
1 GeV. 
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4 Comparison with experiments 



The model was used to perform a fit of the data on structure function F 2 (x, Q 2 ) and 
diffractive structure function (x, Q 2 , in the region of small x (x < 1CT 2 ) and 
Q 2 < 10 GeV 2 [J In Table 1, the full list of parameters (fitted and fixed) is given. 

In Fig. 2, the results for F2(x,Q 2 ) are given as a function of x for different 
values of Q 2 . The description of the data is good. S and L contributions are shown 
separately in order to see the different behaviors. The S contribution is almost 
negligible for very small Q 2 , becoming comparable to the L one at larger Q 2 values. 

In Figs. 3 to 6 we compare the model with the experimental data on diffraction. 
In order to do such a comparison, it is necessary to take into account that different 
experiments use slightly different definitions of diffractive events. In this way we have 
multiplied eq. (Hlf) by a factor Cdiff—1-1 in order to compare with data from 
HI experiment and C(&//=1.3 for ZEUS|. With these factors we take into account 
the different cuts in the mass of the diffractively dissociated proton (larger in the 
case of ZEUS). As in the previous case, we plotted L, S and PPP contributions 

(3) 

separately. In Fig. 3, we show our results for the p— dependence of xpF 2D for 

xp=0.003 and for two values of Q 2 . In Fig. 4, the results are given as a function of 

xp for different values of j3 and Q 2 . For the highest values of Q 2 , only comparison 

with [3=0.4 and /3=0.65 are given. For smaller values of (3, QCD evolution becomes 

important. In Figs. 5, the energy dependence of the diffraction cross-section is 

shown for different values of M and Q 2 . In Fig. 6, the M 2 — dependence of the 

model on diffractive dissociation in photoproduction is compared with HERA data 

for two different energies. Only data with M 2 < 100 GeV 2 are shown for comparison. 
* Actually, for F 2 only data with Q 2 < 3.5 GeV 2 were included in the fit. 

§ Notice that we have taken these values as constant for each experiment, though they could 
also depend on M. This would improve the agreement in Fig. 5.1. 
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For larger values, the effect of the non- diffract ive RRP contribution (not included 
in the model) is expected to be large. 

5 Conclusions 

We have introduced a model of the eikonal type to describe total and diffractive 
7*p interactions. The 7*p interaction is viewed as that of a qq pair, produced by the 
virtual photon, with the proton. The 7*p total cross-section is separated into two 
components : large size (L) for r > r and small size (S) for r < r , where r is the 
transverse distance between q and q. The value of tq - treated as a free parameter - 
turns out to be r ~ 0.2 fm. For the L-component, all the Q 2 -dependence is given 
by the coupling of 7* to the large size qq pair - which is taken as 1/Q 2 at large Q 2 
(eq. (|J)). For the S-component, the Q 2 -dependence is given by the wave function 
of the qq pair (eqs. (0, [3])), computed in perturbative QCD. At large Q 2 , r 2 ~ 1/Q 2 , 
and the unitarity corrections of the S component are higher twist, whereas those of 
the L don't depend on Q 2 . 

A good description of the small x data is obtained both for F 2 and diffractive 
production, in a broad region of Q 2 (0 < Q 2 < 10 GeV 2 ), with a single Pomeron of 
intercept ctp(0) = 1.2. For larger values of Q 2 , QCD evolution becomes important. 
In particular it will give rise to a behavior F<i ~ x~ Ap with Ap significantly larger 
than 0.2 at large Q 2 [13]. For diffraction, this evolution has rather small effects at 



intermediate values of (3 This allows us to use our model, in this case, without 
QCD evolution, up to rather large Q 2 and moderate (3. 

In the region < Q 2 < 10 GeV 2 the unitarity effects are very important and 
produce a significant decrease of the effective Pomeron intercept «p(0) = 1 + Ap 
with decreasing Q 2 . This decrease is controlled by the strength of the unitarity 
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corrections. This, in turn, is controlled by the ratio er'*^" / cr( io *) and its dependence 
on Q 2 . Hence the importance of describing both total cross-sections and diffractive 
production. In our case, Xl > Xs an d the unitarity corrections are more important 
in the L component than in the S. Moreover, these corrections are higher twist at 
large Q 2 in the second case. This is more clearly seen in the diffraction, where the 
S contribution to xpF^(x, Q 2 , ft) is much smaller than the L one for all but the 
larger ft values. 

An important result of our analysis is not only the fact that we can describe 
the data on both structure function and diffractive production in a broad region of 
Q 2 with a single Pomeron, but also that we can describe diffractive production at 
Q 2 = and at intermediate Q 2 using the same value of the triple Pomeron coupling 
(which appears in our parameter a). 

Finally, we would like to discuss the large £n(l/x) limit of the total 7*p cross- 
section in our model. The s,Q 2 ), given by eqs. (pp-fllTD, tend to saturate 
fastly with increasing s to the value 1/2C, due to the large Xl(s, b, Q 2 ). The situa- 
tion is different for the S component. Let's forget for the moment about the triple 
pomeron contribution, i.e. consider the case a=0. As we have said, unitarity correc- 
tions are much smaller in the S component, so, saturation will take place at much 
bigger energies, when the exp(£Ap) term gets large enough. For such energies, cross 
section in the small impact parameter will saturate to a Q 2 -independent value and 
F 2 (x,Q 2 ) ~ Q 2 . This is the usual picture in perturbative QCD 0-0. However, 
by including the nonperturbative (large distance) PPP terms (a ^ 0) we obtain a 
different behavior. Indeed, the large exp(£Ap) factors in the numerator and denom- 
inator of eq. (|T4| ) cancel with each other, and we have a^*p ~ ^2f(£nQ 2 ). Thus, 
the 1/Q 2 smallness of the 7*p cross-section is maintained in the limit x — » 0. 
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Figure captions 



Figure 1. A generic reggeon diagram of our model. It contains the s-channel 
iteration of Pomeron and / exchanges, triple Pomeron (PPP + PfP) diagrams, as 
well as multiple t-channel branchings of the Pomeron of the fan-diagram type. 
Figure 2. F 2 (x,Q 2 ) as a function of x for different values of Q 2 compared with 
experimental data from HI 1995 Jy| (open squares), ZEUS 1995 [jrjj (black circles), 
E665 |17J (black triangles) (notice that the corresponding Q 2 values of these data are 
slightly different) and ZEUS BPT97 |18| (open circles). Dotted curve corresponds 
to the L contribution, dashed one to the S contribution and solid one to the total 
F 2 (x, Q 2 ) given by the model. 

Figure 3. xpF^p as a function of (3 for fixed xp = 0.003 and for Q 2 =4.5 GeV 2 and 
Q 2 =7.5 GeV 2 . Experimental data are from Dotted lines correspond to PPP 
contribution, dashed ones to L term and dotted-dashed to S one. 
Figure 4.1. xpF^p as a function of xp for Q 2 =4.5 and 7.5 GeV 2 and fixed (3=0.04, 
0.1, 0.2, 0.4, 0.65 and 0.9. The curves correspond to the convention of Fig. 3. 



Experimental data are from 119 



Figure 4.2. XpF^o as a function of xp for Q 2 =9, 12 and 18 GeV 2 and fixed (3=0.4 
and 0.65. The curves correspond to the convention of Fig. 3. Experimental data 



are from [19 



Figure 5.1. Energy dependence of diffractive cross section for M=2, 5 and 11 GeV 
and Q 2 =8 and 14 GeV 2 . Experimental data are from [p0[| . The curves correspond 
to the convention of Fig. 3. 

Figure 5.2. Energy dependence of diffractive cross section for different mass in- 
tervals and for low-Q 2 compared with experimental data from [21]. The curves 
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correspond to the convention of Fig. 3. 

Figure 6. Diffractive photoproduction cross-section for W=187 and 231 GeV as 
a function of M 2 from |22| compared with our model. The curves follow the same 
convention as in Fig. 3. 
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Table 1 



Fixed Parameters 


Fitted Parameters 


A P 


0.2 


<?i(o) 


4.56xl0" 3 


A/ 


-0.3 




1.97 GeV" 2 


a' P 


0.25 GeV" 2 




0.56 GeV" 2 


a' f 


0.9 GeV" 2 




0.79 GeV 2 


d2 

^OkL 


3 GeV" 2 


a 


4.63 xlO" 2 GeV" 2 


-"-ops 


2 GeV" 2 


m{ 


0.59 GeV 2 




2.2 GeV" 2 


C s 


0.18 


7/ 


8 


r 


1.06 GeV" 1 


C 


1.5 


m s 


0.15 GeV 2 
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Figure 1 
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Figure 2 
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Figure 3 
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Figure 4.1 
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Figure 4.2 
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Figure 5.1 
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Figure 5.2 
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Figure 6 
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